Relatedness and genomic inbreeding in a sample of Timor ponies.

Timor ponies (TP) were first shipped to Australia in the early 1800s and were highly valued as transport and pack animals, which resulted in TPs contributing to the development of Australian horse breeds. Today, while the exact number of TPs in Australia is currently unknown, there has been recent interest in establishing a domestic breeding program for Australian TPs. The aim of this study was to evaluate the relatedness of a sample of TPs, as well as provide estimates of genomic inbreeding levels to better inform the feasibility of using these animals as founders for a domestic breeding program. Hair samples from each horse were genotyped using the Illumina 80K Infinium Equine genotyping array and data were analysed using PLINK v1.90b7, KING 2.3.2 and R v4.3.1. The results illustrate that there are distantly related and minimally inbred horses within the sampled TPs. Lengths of the ROH segments also indicated that recent inbreeding events are likely to only have occurred in a third of the horses. Overall, these results are promising for the success of a domestic TP breeding program; however, considering the low number of domestic TPs known to reside in Australia, there would certainly still be substantial benefits to incorporating additional TPs either directly from Timor or from areas in Australia that are believed to contain wild descendants of TPs.

Genetic architecture and polygenic risk score prediction of degenerative suspensory ligament desmitis (DSLD) in the Peruvian Horse.

Introduction: Spontaneous rupture of tendons and ligaments is common in several species including humans. In horses, degenerative suspensory ligament desmitis (DSLD) is an important acquired idiopathic disease of a major energy-storing tendon-like structure. DSLD risk is increased in several breeds, including the Peruvian Horse. Affected horses have often been used for breeding before the disease is apparent. Breed predisposition suggests a substantial genetic contribution, but heritability and genetic architecture of DSLD have not been determined. Methods: To identify genomic regions associated with DSLD, we recruited a reference population of 183 Peruvian Horses, phenotyped as DSLD cases or controls, and undertook a genome-wide association study (GWAS), a regional window variance analysis using local genomic partitioning, a signatures of selection (SOS) analysis, and polygenic risk score (PRS) prediction of DSLD risk. We also estimated trait heritability from pedigrees. Results: Heritability was estimated in a population of 1,927 Peruvian horses at 0.22 ± 0.08. After establishing a permutation-based threshold for genome-wide significance, 151 DSLD risk single nucleotide polymorphisms (SNPs) were identified by GWAS. Multiple regions of enriched local heritability were identified across the genome, with strong enrichment signals on chromosomes 1, 2, 6, 10, 13, 16, 18, 22, and the X chromosome. With SOS analysis, there were 66 genes with a selection signature in DSLD cases that was not present in the control group that included the TGFB3 gene. Pathways enriched in DSLD cases included proteoglycan metabolism, extracellular matrix homeostasis, and signal transduction pathways that included the hedgehog signaling pathway. The best PRS predictive performance was obtained when we fitted 1% of top SNPs using a Bayesian Ridge Regression model which achieved the highest mean of R2 on both the probit and logit liability scales, indicating a strong predictive performance. Discussion: We conclude that within-breed GWAS of DSLD in the Peruvian Horse has further confirmed that moderate heritability and a polygenic architecture underlies the trait and identified multiple DSLD SNP associations in novel tendinopathy candidate genes influencing disease risk. Pathways enriched with DSLD risk variants include ones that influence glycosaminoglycan metabolism, extracellular matrix homeostasis, signal transduction pathways.

Population Genetic Analysis of the Estonian Native Horse Suggests Diverse and Distinct Genetics, Ancient Origin and Contribution from Unique Patrilines.

The Estonian Native Horse (ENH) is a medium-size pony found mainly in the western islands of Estonia and is well-adapted to the harsh northern climate and poor pastures. The ancestry of the ENH is debated, including alleged claims about direct descendance from the extinct Tarpan. Here we conducted a detailed analysis of the genetic makeup and relationships of the ENH based on the genotypes of 15 autosomal short tandem repeats (STRs), 18 Y chromosomal single nucleotide polymorphisms (SNPs), mitochondrial D-loop sequence and lateral gait allele in DMRT3. The study encompassed 2890 horses of 61 breeds, including 33 ENHs. We show that the expected and observed genetic diversities of the ENH are among the highest within 52 global breeds, and the highest among 8 related Northern European ponies. The genetically closest breeds to the ENH are the Finn Horse, and the geographically more distant primitive Hucul and Konik. ENH matrilines are diverse and relate to draught and Pontic-Caspian breeds. ENH patrilines relate to draught breeds, and to a unique haplogroup not described before. None of the 33 ENHs carried the "gait" mutation, but the mutation was found in 2 Huculs. The study demonstrates that the ENH is a genetically distinct and diverse breed of ancient origin with no notable pressure of selective breeding.

Genetic diversity and origin of the feral horses in Theodore Roosevelt National Park.

Feral horses in Theodore Roosevelt National Park (TRNP) represent an iconic era of the North Dakota Badlands. Their uncertain history raises management questions regarding origins, genetic diversity, and long-term genetic viability. Hair samples with follicles were collected from 196 horses in the Park and used to sequence the control region of mitochondrial DNA (mtDNA) and to profile 12 autosomal short tandem repeat (STR) markers. Three mtDNA haplotypes found in the TRNP horses belonged to haplogroups L and B. The control region variation was low with haplotype diversity of 0.5271, nucleotide diversity of 0.0077 and mean pairwise difference of 2.93. We sequenced one mitochondrial genome from each haplotype determined by the control region. Two complete mtDNA sequences of haplogroup L were closely related to the mtDNA of American Paint horse. The TRNP haplotype B did not have close matches in GenBank. The phylogenetic test placed this sequence in a group consisting of two horses from China, one from Yakutia, and one from Italy raising a possibility of historical transportation of horses from Siberia and East Asia to North America. Autosomal STR loci were polymorphic and indicated that the TRNP horses were distinctly different from 48 major horse breeds. Heterozygosity, mean number of alleles, and other measures of diversity indicated that TRNP herd diversity was below that observed for most other feral herds and domestic breeds. Both mtDNA and STRs demonstrated that the existing genetic data sets of horses are insufficient to determine the exact origins of the TRNP horses. However, measures of nuclear and mitochondrial diversity have elucidated management needs. It is recommended that new genetic stock be introduced and that adaptive management principles are employed to ensure that unique mitochondrial lineages are preserved and genetic diversity is increased and maintained over time.

The legend of the Canadian horse: genetic diversity and breed origin.

The Canadian breed of horse invokes a fascinating chapter of North American history and as such it is now a heritage breed and the national horse of Canada. The aims of this study were to determine the level of genetic diversity in the Canadian, investigate the possible foundation breeds and the role it had in the development of the US horse breeds, such as Morgan Horse. We tested a total of 981 horses by using 15 microsatellite markers. We found that Canadian horses have high values of genetic diversity indices and show no evidence of a serious loss of genetic diversity and the inbreeding coefficient was not significantly different from zero. Belgian, Percheron, Breton and Dales Pony, unlike the light French horses, may have common ancestries with the Canadian and could be important founders. However, the Shire and Clydesdale influenced the Canadian to a lesser extent than French and Belgian draft breeds. Furthermore, our finding indicated that there was no evidence of a clear relationship between Canadian and Oriental or Iberian breeds. Also, the Canadian likely contributed to the early development of the Morgan. Finally, these findings support the ancient legends of the Canadian Horse as North America's first equine breed and the foundation bloodstock to many American breeds and may help in the management and breeding program of this outstanding breed in North America.

Copy number variation in the horse genome.

We constructed a 400K WG tiling oligoarray for the horse and applied it for the discovery of copy number variations (CNVs) in 38 normal horses of 16 diverse breeds, and the Przewalski horse. Probes on the array represented 18,763 autosomal and X-linked genes, and intergenic, sub-telomeric and chrY sequences. We identified 258 CNV regions (CNVRs) across all autosomes, chrX and chrUn, but not in chrY. CNVs comprised 1.3% of the horse genome with chr12 being most enriched. American Miniature horses had the highest and American Quarter Horses the lowest number of CNVs in relation to Thoroughbred reference. The Przewalski horse was similar to native ponies and draft breeds. The majority of CNVRs involved genes, while 20% were located in intergenic regions. Similar to previous studies in horses and other mammals, molecular functions of CNV-associated genes were predominantly in sensory perception, immunity and reproduction. The findings were integrated with previous studies to generate a composite genome-wide dataset of 1476 CNVRs. Of these, 301 CNVRs were shared between studies, while 1174 were novel and require further validation. Integrated data revealed that to date, 41 out of over 400 breeds of the domestic horse have been analyzed for CNVs, of which 11 new breeds were added in this study. Finally, the composite CNV dataset was applied in a pilot study for the discovery of CNVs in 6 horses with XY disorders of sexual development. A homozygous deletion involving AKR1C gene cluster in chr29 in two affected horses was considered possibly causative because of the known role of AKR1C genes in testicular androgen synthesis and sexual development. While the findings improve and integrate the knowledge of CNVs in horses, they also show that for effective discovery of variants of biomedical importance, more breeds and individuals need to be analyzed using comparable methodological approaches.

Genetic diversity in the modern horse illustrated from genome-wide SNP data.

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Genome-wide analysis reveals selection for important traits in domestic horse breeds.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.

Maternal lineages in native Canadian equine populations and their relationship to the Nordic and Mountain and Moorland pony breeds.

A 378-bp section of the mitochondrial displacement loop was used to estimate genetic diversity in the native Canadian equine populations. The inclusion of 10 Mountain and Moorland, 3 Nordic pony breeds, 2 feral populations, and 5 horse breeds were also investigated as they may have influenced the development (or rejuvenation) of the native Canadian populations. A total of 281 samples were sequenced, which produced 75 haplotypes derived from 54 informative sites. On further investigation, 36 of these 75 haplotypes were found to be previously unreported. Overall, total diversity was lowest in the feral Sable Island population with a haplotype diversity (0.27 ± 0.12), nucleotide diversity (0.0007 ± 0.0004), and pairwise difference of 0.286 ± 0.317. This is not surprising due to the geographic isolation of this population. Haplotype diversity was highest (1.00 ± 0.13) in the New Forest population, pairwise difference was highest (8.061 ± 4.028) in the Icelandic breed, whereas nucleotide diversity was highest in the Exmoor breed (0.0209 ± 0.0025). Within the Canadian populations, haplotype diversity was highest in the Newfoundland pony (0.96 ± 0.08), whereas pairwise difference and nucleotide diversity was highest in the Canadian horse (7.090 ± 3.581 and 0.0188 ± 0.0042, respectively). Three different estimates of genetic distances were used to examine the phylogenetic relationships amongst these populations. All 3 estimates produced similar topologies. In general, the native Canadian populations were highly represented in the D clade, with particular emphasis in the D1 and D2 clades. This is an important factor when considering the phylogenetic conservation of these Canadian equine populations.

Genetic variability in the Skyros pony and its relationship with other Greek and foreign horse breeds.

In Greece, seven native horse breeds have been identified so far. Among these, the Skyros pony is outstanding through having a distinct phenotype. In the present study, the aim was to assess genetic diversity in this breed, by using different types of genetic loci and available genealogical information. Its relationships with the other Greek, as well as foreign, domestic breeds were also investigated. Through microsatellite and pedigree analysis it appeared that the Skyros presented a similar level of genetic diversity to the other European breeds. Nevertheless, comparisons between DNA-based and pedigree-based results revealed that a loss of genetic diversity had probably already occurred before the beginning of breed registration. Tests indicated the possible existence of a recent bottleneck in two of the three main herds of Skyros pony. Nonetheless, relatively high levels of heterozygosity and Polymorphism Information Content indicated sufficient residual genetic variability, probably useful in planning future strategies for breed conservation. Three other Greek breeds were also analyzed. A comparison of these with domestic breeds elsewhere, revealed the closest relationships to be with the Middle Eastern types, whereas the Skyros itself remained isolated, without any close relationship, whatsoever.

Genetic variability in the Skyros pony and its relationship with other Greek and foreign horse breeds

In Greece, seven native horse breeds have been identified so far. Among these, the Skyros pony is outstanding through having a distinct phenotype. In the present study, the aim was to assess genetic diversity in this breed, by using different types of genetic loci and available genealogical information. Its relationships with the other Greek, as well as foreign, domestic breeds were also investigated. Through microsatellite and pedigree analysis it appeared that the Skyros presented a similar level of genetic diversity to the other European breeds. Nevertheless, comparisons between DNA-based and pedigree-based results revealed that a loss of genetic diversity had probably already occurred before the beginning of breed registration. Tests indicated the possible existence of a recent bottleneck in two of the three main herds of Skyros pony. Nonetheless, relatively high levels of heterozygosity and Polymorphism Information Content indicated sufficient residual genetic variability, probably useful in planning future strategies for breed conservation. Three other Greek breeds were also analyzed. A comparison of these with domestic breeds elsewhere, revealed the closest relationships to be with the Middle Eastern types, whereas the Skyros itself remained isolated, without any close relationship, whatsoever.

Inbreeding and genetic structure in the endangered Sorraia horse breed: implications for its conservation and management.

The Sorraia horse is a closed breed with reduced effective population size and considered in critical maintained risk status. The breed exists in 2 main breeding populations, one in Portugal and one in Germany, with a smaller population size. A set of 22 microsatellite loci was used to examine genetic diversity and structure of the Sorraia horse breed and to compare individual inbreeding coefficient F, estimated from pedigree data, with individual heterozygosity and mean d(2). The Sorraia horse shows lower levels of microsatellite diversity when compared with other horse breeds. Due to management strategies, there are clear differences in the genetic structure of the 2 main Sorraia horse populations. Individual heterozygosity was shown to be a good estimator, used together with or as an alternative to inbreeding coefficient, in predicting fitness and evaluating the inbreeding level of the Sorraia horse. The information gathered in this study, combined with information available from previous studies, offers an important and wide information base for the future development of an effective breeding management of the Sorraia horse in order to preserve this endangered breed.

Genetic structure and phylogenetic relationships of the Polish Heavy horse.

In this study a wide range of genetic markers (12 microsatellites, 7 blood-group loci, 10 blood-protein loci) and mitochondrial DNA (mtDNA) were used to assess genetic diversity in Polish Heavy horses. Three random samples were sequenced for 421 bp of the mitochondrial D-loop region, but no clear phylogenetic patterns were seen in mtDNA variation. Both heterozygosity and diversity levels are fairly high in Polish Heavy horses. In phylogenetic analysis the draught horses form a distinct cluster that pairs with the true pony breeds. Within this 'cold-blooded' group, the Polish Heavy Horse clusters most closely with the Posavina breed from Croatia and the Breton breed from France. From the standpoint of genetic conservation, the Polish Heavy Horse does not appear to be in jeopardy.

Iberian origins of New World horse breeds.

Fossil records, archaeological proofs, and historical documents report that horses persisted continuously in the Iberian Peninsula since the Pleistocene and were taken to the American continent (New World) in the 15th century. To investigate the variation within the mitochondrial DNA (mtDNA) control region of Iberian and New World horse breeds, to analyze their relationships, and to test the historical origin of New World horses, a total of 153 samples, representing 30 Iberian and New World breeds, were analyzed by sequencing mtDNA control region fragments. Fifty-four haplotypes were found and assigned to seven haplogroups. Reduced levels of variation found for the Menorquina, Sorraia, and Sulphur Mustang breeds are consistent with experienced bottlenecks or limited number of founders. For all diversity indices, Iberian breeds showed higher diversity values than South American and North American breeds. Although, the results show that the Iberian and New World breeds stem from multiple origins, we present a set of genetic data revealing a high frequency of Iberian haplotypes in New World breeds, which is consistent with historical documentation.

Mitochondrial DNA D-loop sequence variation among 5 maternal lines of the Zemaitukai horse breed

Genetic variation in Zemaitukai horses was investigated using mitochondrial DNA (mtDNA) sequencing. The study was performed on 421 bp of the mitochondrial DNA control region, which is known to be more variable than other sections of the mitochondrial genome. Samples from each of the remaining maternal family lines of Zemaitukai horses and three random samples for other Lithuanian (Lithuanian Heavy Draught, Zemaitukai large type) and ten European horse breeds were sequenced. Five distinct haplotypes were obtained for the five Zemaitukai maternal families supporting the pedigree data. The minimal difference between two different sequence haplotypes was 6 and the maximal 11 nucleotides in Zemaitukai horse breed. A total of 20 nucleotide differences compared to the reference sequence were found in Lithuanian horse breeds. Genetic cluster analysis did not shown any clear pattern of relationship among breeds of different type

Variation in the mitochondrial control region sequence between the two maternal lines of the Sorraia horse breed

The mitochondrial control region sequence was compared between the two surviving matrilineal lineages of the Sorraia horse. Two haplotypes were found, showing three differences between them. The reduced number of surviving maternal lineages emphasizes the importance of establishing a conservation plan for this endangered breed